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ABSTRACT 

The expansion of the universe may be observed in "realtime" by measuring changes in the patterns of the 
anisotropy in the CMB. As the universe ages, the surface of decoupling — or the CMB photosphere — moves 
away from us and samples a different gravitational landscape. The response of the CMB to this new landscape 
results in a different pattern than we observe today. The largest change occurs at £ « 900. We show that with an 
array of detectors that we may envision having in a couple of decades, one can in principle measure the change 
in the anisotropy with two high precision measurements separated by a century. 
Subject headings: cosmic microwave background - cosmology: observations 



1. INTRODUCTION 

Measurements of the anisotropy of the cosmic microwave 
background (CMB) provide a snapshot of the universe some 
380,000 years after the big bang, when the primordial plasma 
decouples from the baryons. The power spectrum of the 
anisotropy is computed with a line of sight integration of the 
coupled Boltzmann equati ons between now and s ome time 
well before decoupling (Seliak & Zaldarriaga 19961) (here- 
after paper 1). This formalism may be extended to consider 
future times when the decoupling surface is further away from 
us. We may think of this as investigating an expanding CMB 
photosphere. Of course we do not know the details of the 
gravitational landscape beyond the current photosphere and 
so we can only probe the future in a statistical sense. 

The conformal time of the decoupling surface in a flat uni- 
verse is obtained from the Friedmann equation as follows: 



r(fl) = 



da' 



(1) 



where a is the scale factor, //q is the current Hubble parameter, 
J7r,o is the current radiation density relative to the critical den- 
sity, ilm.o is the current matter density relative to the critical 
density, and f^A.o is the current dark energy density relative to 
the critical density. At all times ilr+ft„,+^A = 1- With param- 
eters inspir ed by the Wilkinson Microwave Anisotropy Probe 
(WMAP) ( Spergeletal. 20071) . we take n,. = 8.056 x 10-\ 
n,„ = 0.257, Oa = 0.743 and Hq = 72 km/s/Mpc. The cur- 
rent conformal distance to the decoupling surface (a = 1 in 
the above equation) is cr = 14373 Mpc. 

2. COMPUTING THE POWER SPECTRUM 

To perform our study, we use CAMB dLewis et al. 20001) 
which follows the algorithm in Selj ak and Zaldarriga's CMB- 
FAST (ISeliak & Zaldarriaga 19961) . We consider only scalar 
perturbations. We modified the CAMB code so that it com- 
putes the following evolution function (eq. 13 in paper 1) 



ATi(Tf,k)= S(T,k)je[k(Tf-T)]dT 



(2) 



where S is the source function, ji is a spherical Bessel func- 
tion, T is the conformal time, and k is the comoving wavevec- 
tor of the perturbation. The modification allows for arbitrary 
final T (as opposed to fixing ry = tq). The source function is 



Table 1. Time slices and Scale Factors 



Time slice 


CT 


Scale factor 


Years 


Mpc 


a 





14373 


1.0000 


1 X 10* 


14403 


1.0074 


5 X 10* 


14523 


1.0372 


1 X 10' 


14669 


1.0754 


2 X lO' 


14944 


1.1545 


4x 10' 


15440 


1.3255 


5 X 10' 


15663 


1.4182 


6x 10' 


15872 


1.5163 


8 X 10' 


16251 


1.7304 


1 X 10'" 


16583 


1.9714 


4x 10'" 


18657 


13.340 


1 X 10" 


19011 


601.95 



Note. — To find the physical time of any time slice add 13.66 X lO' 
years, the current age of the universe for our adopted parameters. Because 
dr = dt/a, the physical time is computed from equation [T] except with an 
additional factor of I /a' in the integrand. Note that the universe doubles in 
diameter when the physical time is 75% greater than today and the conformal 
time is 15% greater 



given by 



S(k,T) = e 



Kf 



d 



g (At-o + + — g-^^ + U>+ * +P(k, T) 



dr k 

(3) 

where P(k,T) represents polarization terms (see eq. 12 in pa- 
per 1) which we ignore for simplicity. Here g(r) = -/ie"" 
is the visibility function with k{t) the optical depth from tq 
(note that for r/ > tq, k/ < 0). Eq. (|3]l differs from eq. 12 
in paper 1 only in the overall e'^f term which accounts for 
the changing optical depth. The first order term in the ex- 
pansion of the Boltzmann equation, A70 (eq- 3a in paper 1), 
the potentials ^I^ and $, the velocity of the baryons vi,, and the 
polarization terms are computed by CAMB and need no mod- 
ification. The primary modification to CAMB is to calculate 
S out to Tf, rather than to tq.' In addition, some of CAMB's 
optimizations that led to sparse sampling at recent times were 
removed. 

' We accomplish this simply by setting CAMB's internal tq variable to the 
maximum Tf in which we are interested. This is valid because in the code the 
concept of the "present" is linked to a = 1 (which we do not modify), rather 
than to T = TQ. 
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To compute the power spectrum at any time in the future 
we simply form the analog of eq 9 in paper 1 : 

C,{Tf) = (47rfJ k^P^(k) [AT,(0,rf,k)A*JO,Tf,k)] dk (4) 

where P\!j(k) is the initial power spectrum. Figure [1] shows 
the power spectrum for future times. We see three major ef- 
fects: (1) The power spectrum amplitude drops off due to the 
1/a scaling of the CMB temperature;^ (2) The features shift 
to smaller angular scales due to the recession of the surface 
of last scattering; (3) The low-i tail becomes enhanced com- 
pared to the peak due to the integrated Sachs-Wolfe (ISW) 
effect caused by the shift to dark-energy dominance. The en- 
hanced ISW effect is not present in runs of the code without 
dark energy. In roughly 25 x 10^ years, the ISW effect at £ = 2 
will exceed the height of the first acoustic peak. 

6000 I 




1 10 100 1000 

Multipole 

Fig. 1 . — The temperature angular power spectrum of the CMB at several 
representative time steps into the future. The scale factors and conformal 
distances are given in Table 1 . The shift to the right of the first acoustic peak 
for some future time is if = ^oC'''/' ~'^rffr)/('n) "^rfcc) where tj^^. = 286 Mpc is 
the conformal time of decoupling, and £o = 220 the current £ of the acoustic 
peak. 



3. MAPS OF THE FUTURE SKY 

To visualize the evolution of the CMB photosphere, we 
make maps of the sky for various physical time slices (see Ta- 
ble 1). For any particular map, which we take as a Gaussian 
random field, all the information is contained in the power 
spectrum, Q =< fl£/„a|,„ >, where the ai„s are the coefficients 
of a spherical harmonic decomposition of the map. To gener- 
ate a visualization, one draws £+1 random normal complex 
deviates, r, with variance Q to generate a set of 2^+ 1 af„,s 
that satisfy ai-m = a}^. From these one forms the real valued 

T(e, (/)) = J2ejn aimYim(0, </>)■ 

A set of maps that shows the future evolution of the CMB 
will be correlated. To account for this, we compute the full 
covariance matrix 

Cy = (47TfJ k^P^(k) [Aj-,(0,T'>)A*^(0,r>,^)] dk, (5) 

where t' is the conformal time at any time in the future. Thus 
for a sequence of say n maps, we would compute anx n ma- 
trix for each £. 

^ Note that Eqs. j2) and (4) are dimensionless (since At is dimension- 
less by definition). We give units to Ci{Tf) by multiplying by /cP'iTf), 
accounting for this effect. 



We now extend the method given above to generate a set of 
correlated fl/,„s. The first step is to decompose the covariance 
matrix as Ci = where is a diagonal matrix. We 

then compute a^,„ = M£[-\/nf ]r where r is a vector of complex 
random deviates. This a£„, has the covariance matrix in eq.|5] 
Figures |2]|5] show a set of four time slices starting with a ran- 
dom full-sky map that follows the WMAP parameters.^ One 
can see that most of the change occurs on small angular scales 
where the photosphere more quickly samples different poten- 
tial wells as it expands. One also sees that as time progresses 
large angular scale fluctuations become more prominent as 
dark energy dominates the expansion. 




-450 ^^^^^^^^^^^-^ 450 



Fig. 2. — Present map for a set of random variables that follows the WMAP 
power spectrum. The units on the color bar are /iK. 




-419 ^^^^^^^^^^^-^ 419 



Fig. 3. — The same map as in Figure[2]but now for a time slice 1 X lO' 
years in the future. Note that the color bar has been rescaled. The primary 
difference between the two maps is at small angular scales. 




-317 ^ 317 



Fig. 4. — The same map as in Figure[2]but now for a time slice 5 X lO' 
years in the future. Again, the color bar has been rescaled. 

Figure|6]shows elements of the covariance matrix as a func- 
tion of time. As expected, nearby time slices are strongly cor- 
related. As time progresses the covariance between current 
and future time slices at small angular scales disappears first, 

' One could start with the WMAP sky though we have not done this. 



Expanding CMB Photosphere 
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Fig. 5. — The same map as in Figure[2]but now for a time slice 40 X lO' 
years in tlie future. Tlie color bar has been rescaled. Note that the large 
angular scale fluctuations are more prominent relative to the smaller angular 
scales. 
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and then later the covariance at large angles decreases. The 
time it takes for the scale factor to double, 10^°y, gives a 
characteristic time for the future sky to decorrelate with what 
we observe today. 

Figure|7]gives the correlation of several time steps with the 
present. In this plot, the decrease in temperature has been 
scaled out so that features in the anisotropy may be compared 
directly. At high £, features in the sky at late times are uncor- 
rected with those at present. However, at low £, the late-time 
ISW features remain correlated for > 100 billion years, indi- 
cating that these are very long lived structures. 
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Fig. 6. — The C^^ element of the covariance matrix for several time slices. 
The power spectrum at to is also plotted for comparison. Some of the de- 
crease in the covariance is attributable to the decrease in the CMB tempera- 
ture. For example, in lO' years the power at £ = 200 decreases to 36% the 
current value but 16% of the decrease arrises because T = T(i/a. 



4. MEASURING THE CHANGE IN THE CMB 

Measuring the difference between two high precision maps 
of the anisotropy taken a century apart offers, in principle, 
a way to directly observe the expansion of the universe. Un- 
like a measurement of the temperature of the CMB, the differ- 
ence between two maps is moderately insensitive to calibra- 
tion. Rather, it is a change in spatial structure that is observed. 
Thus one needs a well understood pointing solution which is 
technically straightforward to achieve. 

Figure [8] shows the power spectrum of the difference be- 
tween two maps taken 100 years apart. We use the formal- 
ism in Knox (1995) to compute the experimental uncertainty. 
Since there is just one sky realization, cosmic variance is ig- 



FlG. 7. — The correlation of several time slices with the present. 



nored and the uncertainty per £ is 



^'^^ = \ f ,^,^,, <^'»/>-.vexp(^^g?/2/81n2), (6) 

y Jsky(^(-+ I) 

where fsky is the fraction of sky covered, a is the uncertainty 
per sky pixel of solid angle Vlpi^, and 0i/2 is the FWHM of the 
beam profile. For the uncertainty bands shown in Figure [8] 
we imagine an array of 3000x3000 detectors at 150 GHz each 
with a sensitivity of 40 mKs'/^ jBock et al. 2 006) observing 
the sky with 0.86' angular resolution. The observations would 
last 4 years, cover 75% of the sky, and would have to be done 
from a satellite. The only element not aheady demonstrated 
is large array. Currently, arrays of 32x32 detectors are being 
built. 




Multipole (I) 

Fig. 8. — The top curves shows the current CMB EE polarization spectrum 
for the parameters in this paper and a cuirent optical depth of 0. 1 . The two 
curves below that show the BB spectrum for tensor to scalar ratios of r = 0.3 
and r = 0.01. The sohd curve at the bottom is the power spectrum of the 
difference between two CMB anisotropy maps made 100 years apart. The 
dotted line is Equation 6 and the hashed boxes show the errors after binning 
in I. 

The fundamental limit to such a measurement is likely to 
be variable point sources and variable foreground emission. 
Though in principle these can be identified and removed spec- 
trally, this capability would add complexity to the "simple" 
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scheme outlined above. For comparison, there are currently 
experiments being designed with the sensitivity to measure 
the r = 0.01 B-mode polarization. The improvement to go 
from these planned missions to measuring the signal we de- 
scribe is on the order of the improvement between observa- 
tions of the 1980s and the current observations. 



We gratefully acknowledge discussions with Mike Nolta, 
Uros Seljak and Paul Steinhardt. Uros provided a key in- 
sight for computing and manipulating the covariance ma- 



trix in equation |5] In the course of this work, we learned 
that Adam Moss, Douglas Scott, and Jim Zibin were work- 
ing on a more extensive version of a related calcula- 
tion. A movie of the expanding photoshpere is available at 
http://phy-page-g5.princeton.edu/^page This paper is based 
on Stuart Lange's senior thesis (Lang^2007) and made exten- 
sive use of the publicly available HEALPix and CAMB soft- 
ware packages. The research was supported by NASA award 
LTSA03-000-0090 and by NSF Grant No. 0355328. 
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